Universality of the momentum band density of periodic networks 
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The momentum spectrum of a periodic network (quantum graph) has a band-gap structure. We 
investigate the relative density of the bands or, equivalently, the probability that a randomly chosen 
momentum belongs to the spectrum of the periodic network. We show that this probability exhibits 
universal properties. More precisely, the probability to be in the spectrum does not depend on 
the edge lengths (as long as they are generic) and is also invariant within some classes of graph 
topologies. 
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The spectrum of Schrodinger operator in periodic 
medium is calculated using the Floquet-Bloch procedure 
PQ: the periodic medium is replaced with its fundamen- 
tal domain endowed with parameter-dependent quasi- 
periodic boundary conditions. The resulting parameter- 
dependent spectrum is called the dispersion relation, and 
the range of the dispersion relation is precisely the spec- 
trum of the original structure. The spectrum has a band- 
gap structure and knowing the band location and sizes is 
of utmost importance in the theories of condensed matter 
and of dielectric and acoustic media [2-6 . Of particular- 
recent interest is understanding the spectrum of quantum 
graphs E] , motivated by their application to solid state 
[9] HO], photonic crystals [11], carbon nano-structures [12] 
as well as their use as models for quantum chaos, both 
in theoretical |13ffT5] and experimental [H?Il2T)] studies. 

In the present Letter we explore the relative size of 
bands and gaps and discover a curious universality. To 
be more precise, we ask the following question: what is 
the probability p s that a randomly and uniformly chosen 
momentum belongs to the spectrum of the graph? For 
example, consider the Z^periodic graphs of Fig. [l] How 
does p a change if we change the lengths in the fundamen- 
tal cell of the graph, from Fig. [TTb) to Fig. Hie)? How 
does p a change if we change the topological structure to 
Fig. gd) or ge)? 

Denote by p a (K) the probability of a uniformly cho- 
sen momentum k G [0, K] to be in the spectrum and let 
p a := lim^^oo p,j (K). We find that the probability p a is 
well-defined and is independent of many features of the 
fundamental cell. In particular, all choices in Fig. flTb) 
to (d) lead to the same value of p a (assuming a generic 
choice of edge lengths). This is illustrated by a numerical 
simulation in Fig. [2] We will derive the limiting value 
analytically below. Note that the value of p a for the cell 
in Fig. [ITe) turns out to be different from the others and 
will also be calculated. 

Let us put the discussion onto a more formal footing. 
We consider a Z d -periodic network of quantum wires on 
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Figure 1: (a) An example of a Z -periodic graph and (b) its 
fundamental cell; (c)-(e) are other examples of the fundamen- 
tal cell. 
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Figure 2: (color online) Numerical simulation of the conver- 
gence of p a (K) is shown for the three periodic graphs from 
Fig. l(b)-(d). The limiting value, p a , is shown as a dashed 
line. The graph lengths are normalized such that K equals 
the average number of spectral bands. 



which we are solving the spectral problem 
d 2 tp 



dx 2 



k 2 vj, (1) 

subject to the Kirchhoff-Ncumann vertex conditions 
ip{x) is continuous at v, 



(2) 



where the sum is over the edges £ v emanating from the 
vertex v and the derivatives are taken into the edge. We 
denote by a the set of k values for which there is a so- 
lution to (flj-pl; this is the momentum spectrum of the 
graph. Now, the definition of p a can be formally written 
as 



Pa = lim p a (K) = lim — 

A'— s-oo K— ¥00 A 



a- n [o, K] 



(3) 



In this Letter we establish several properties of the 
probability p a - First of all, the above limit always exists. 
In addition, if there is at least one gap in the spectrum, 
there are infinitely many gaps and p a < 1. Similarly, if 
there is at least one non-flat band, there are infinitely 
many and p a > 0. Finally, and perhaps most strikingly, 
provided the lengths of edges in the fundamental set are 
generic, the value of p a is independent of their precise 
value. We also find that the value of p a is independent 
of some details of the cell's topology. 

Secular equation and dispersion relation. In the 
Floquet-Bloch procedure for quantum graphs (see, e.g., 
[7]) we identify a set of J generators of the lattice of pe- 
riods and assign to each a quasi-momentum variable <x,-, 
j = 1, . . . , J. If the vertices i>+ and w_ of the fundamen- 
tal cell are identified by the action of the j'-th generator, 
we impose the quasi-periodic conditions 

t/j( v+ ) = e ia ^(v-), i>'(v+) = -e ia *<ip'(v-). (4) 

We remind the reader that we use the convention of al- 
ways taking the derivatives into the edge, which explains 
the minus sign in conditions (El. For example, in the fun- 
damental cell of Fig. flTb) the empty circles denote the 
vertices connected through the condition of the above 
type. Identifying these periodically related vertices cre- 
ates new cycles, Cj, j = 1, . . . , J, on the graph and the 
resulting problem is equivalent to a graph with magnetic 
fluxes aj through the corresponding cycles. For exam- 
ple, the result of the Floquet-Bloch procedure for the 
fundamental cell in Fig. fl[d) is equivalent to the mag- 
netic graph in Fig. pla). We denote by E the number of 
edges of the resulting magnetic graph. 

Expanding the solutions to (fTl) in the basis of e ±tkx and 
applying the vertex conditions leads, after some linear 
algebra (see [13]), to the secular equation 



F(k;a) :=det(l-e i(A+feL >S 



0, 



(5) 



where all matrices act in the space of coefficients on di- 
rected edges; each edge gives rise to two directed edges of 
equal length, therefore all matrices have degree 2E. The 
diagonal matrix L is the matrix of lengths of the directed 
edges. The diagonal matrix A contains the magnetic 
fluxes ctj that are put upon the edges created by vertex 
identifications. The magnetic fluxes change sign when re- 
versing the direction of the corresponding edge. Finally, 





(a) 



(b) 



Figure 3: (a) A graph consisting of a loop pierced by a mag- 
netic flux and a decoration (b) Similar graph, but with a single 
edge decoration 



the unitary matrix S contains directed edge-to-edge scat- 
tering coefficients, which, for scattering at a Neumann- 
Kirchhoff vertex of degree d, is equal to — 1 + 2/d for 
back-scattering and 2/d for forward scattering. Most im- 
portantly, for our vertex conditions the matrix S is inde- 
pendent of k. See (10 1— (111, which show these matrices 
for a specific graph. 



Next we apply a clever trick originally due to Barra and 
Gaspard [2T] (see also |Hj): we introduce a new function 
<£>(«; a) such that 



$(ki = kli, ■ . . , ke — Me; a) := F(k; a), 



(6) 



where Ii,...,Ie are the graph edge lengths. A cursory 
look at equation (pf reveals that the variables n e , e = 
1, . . . , E need only be known modulo 2ir. For a fixed a, 
define XU to be the set of solutions of 



<&(k;o?) = 0, 



(7) 



on the torus T^- := [0, 2ir) . Then the roots k n of the 
equation F(k; a) = can be interpreted as the times (k 
values) of piercing of the set Eg by the flow 



K (k) = k ■ (hjfo, ■ ■ • , Ie) mod27r. 



(8) 



We now conclude that k belongs to the spectrum, a, of 
the periodic graph if the corresponding point n(k) be- 
longs to the set £5 for some value of a (which itself 
belongs to a J-dimensional torus). For future purposes 
we define 



U £* 



(9) 



5<E[0,2tt) 



We will now compute the set E in a simple but important 
example and then proceed to discuss how the questions 
about the band probability p a can be related to the prop- 
erties of the set E. 

Loop with an edge. We now compute the set E for a 
graph which consists of a loop pierced by magnetic field 
with flux a and a single edge attached, see Fig. I3jb) . 

The numbering of the directed edges is given in 
Fig. p[b). According to this numbering the matrices A, 
L and S are given by 



A = diag(a, 0, -a, 0), L = diag(Zi, h, h, h 



(10) 




Figure 4: (color online) (a) The zero sets, E a , of $ («i, K2; a) 
are shown for a range of values a £ [0, it] using a blue (a — 0) 
- red (a — n) color scale, (b) The set E = Uaefo 2ir) ^ a 
is shaded, its blue boundaries are the zeros of <E> (•, •; a — 0) 
and the red boundaries are the zeros of $ (•, •; a = 7r). A flow 
k (k) = k ■ (l\, I2) on the torus is indicated. The bands of the 
spectrum a are the solid black segments of the flow line; the 
gaps are drawn in light gray. 



and 



' 2/3 2 /3 - 1 / 3 0" 

1 

-1/3 2 /3 2 /3 

. 2 /3 -V3 2 /3 0, 



(11) 



The secular function $ evaluates to (up to some non- 
zero factors) 



$ = 2cos(k 2 ) (cos(ki) 



cos(a)) 

— sin(Ki) sin(K 2 ). (12) 



The zero sets E a for a range of values of the parameter a 
are shown on Fig. pja) . Note that it is enough to consider 
the values a € [0, w] as E_ Q = E Q (see (12l). 

Probability to be in the spectrum. From the discussion 
above we conclude that the probability p a for a random k 
to be in the spectrum a is equal to the proportion of time 
the flow defined by (8| spends in the set E. Depending 
on the commensurability properties of the set of the edge 
lengths, {l e } e= i, the flow covers densely the entire torus 
or is restricted to a flat submanifold 



L := span (A; • (h,l2, ■ ■■) modTE : k G R) 

= {xGlR B : Mx = 0} modT B , 

where M is a matrix with rational coefficients (it gives the 
rational dependencies in the length sequence (Zj, . . . , Ie))- 
In the latter case, the flow is ergodic on the submanifold 
L. The probability p a is therefore the relative volume 



Pa 



vol L (L n E) 
vol L (i) : 



(13) 



where the subscript L indicates that the volume should 
be taken in the appropriate dimension (equal to E minus 
the rank of the matrix M). Formula (13 1 remains valid 



in the case of rationally independent lengths, when we 
simply take L to be the entire torus. This immediately 
implies that the probability p a remains the same as long 
as the edge lengths are rationally independent. 

Returning to our example, we calculate p a explicitly. 
Using symmetry we compute the area of 1/8-th of the 
set E, the part in the lower left corner. It is bounded by 
the coordinate axes and the set E w , which from (12) we 
re-parameterize as 



tan(«2) = 2cot(«i/2). 



Therefore the ratio in ( 13 1 evaluates to 



Pa 



tan -1 (2 cot («/a)) dn sa 0.64. 



(14) 



(15) 



We can further prove that this universality of p a ex- 
tends to a certain class of decoration structures. These 
are the decorations that attach to the base line by means 
of a single edge, as in Fig.JTJa) to (d). Proving the univer- 
sality is done by reducing the influence of the decoration 
on the secular equation to a single scattering reflection 
phase located at the degree one vertex of the graph in 
figure [3Tb). The phase enter the matrix S as follows, 



S = 



2 /3 2 / 3 - 1 / 3 

3(k 3 ,...,k e ) 

-1/3 2 /3 2/3 

2/ 3 _l/ 3 2/3 



(16) 



While the precise form of the phase («3, . . . , Ke) may 
be complicated, its effect on the function <& gets averaged 
out by ergodicity. More precisely, we now assume that 
the rational relations (if any) defining the submanifold L 
do not involve k,\ and k.2- In other words, the lengths 
of the edges 1 and 2 are rationally independent of each 
other and of the lengths of the decoration's edges. We 
need not assume anything about the lengths of edges of 
the decoration. 

One can now easily read from the de- 
terminant (see ((sl,©) that the function 
$ has the form $ (m, K2, ■ . ■ , Ke'> cc) = 



$(ki,k 2 + 1/26(^3, 



the RHS is as in (12 1. 
variables 



,ke)',ch), where $(•,•; a) in 
Introducing the change of 



k 2 



K 2 + -6(k 3 ,...,ke) 



(17) 



the integrals in (13 1 factorizc. Namely, denote by T 2 the 
torus with respect to k\ and £2 and by Te-2 the torus 
with respect to the other variables. Note that the set E 
depends only on the variables Ki and £2 (and is cylindri- 
cal with respect to the other variables) . The submanifold 
L, on the other hand, is cylindrical with respect to K\ and 
hi- Therefore 



Pa 



vo1t 2 (Tg n E) vo1 Te _ 2 (T £ _ 2 n L) 
vo1t 2 (T 2 ) vo1 Te _ 2 (Je-2 n L) 



(18) 



reducing to the expression in ( 13 1 , where L there is identi 



fied as T2 in ( 18 ). We thus proved that for all decorations 



of the type discussed above the probability to be in the 



spectrum is given by (151 



To give a final example of a different nature, for the 
fundamental cell depicted in Fig. [IJe), the secular equa- 
tion can be shown to be equivalent to 

, 1 . . . 

sm(K! + K2 + K3) sin K\ sin k 2 sm k 3 

= sinKi + cosa(sinK 2 + sinK 3 ), (19) 

and the corresponding value of p a was calculated numer- 
ically to be 0.43. 

Conclusions. The arguments presented above apply 
to all graphs and result in three general conclusions. 
First, given a Z d -periodic graph with an arbitrary fun- 
damental cell, the probability p a is independent of the 
specific edge lengths, as long as there are no rational de- 
pendencies between some of them. Even if such depen- 
dencies exist, an appropriate ergodicity argument shows 
that the limit ((3| which defines p a exists and its exact 
value depends on the nature of the edge lengths rational 
dependencies (as well as the graph's topology) . Secondly, 
we have shown that p a is robust even within some topo- 
logical modifications of the graph - attaching a prescribed 
class of decorations. Thirdly, if there exists at least one 
non-flat band (resp. gap) in the spectrum, it must arise 
from an open set on the torus which is a subset of £ 
(corresp. T\S). The ergodic flow on the torus will pass 
through this set infinitely many times, resulting in an infi- 
nite number of non-flat bands (resp. gaps) of comparable 
size. From equation ( 13 1 we can immediately conclude 
that p a > (resp. p a < 1). 

Our setup calls for comparison with periodic poten- 
tials on the line, in particular the singular potentials 8 
and 8' [23]. Note that we measure our band and gap 
sizes in terms of the momentum variable k, not energy 
(which scales as k 2 ). For smooth periodic potentials and 
8 potentials, the gaps sizes decrease as k — > 00, while 
the band lengths converge to a constant, resulting in 
Pa = 1 [24"ll25| . The 6' potential has an opposite behav- 
ior, asymptotically equivalent to disconnecting the graph: 
the band lengths decrease and the gaps approach a con- 
stant size, resulting in p ff =0 [25 . Our results show 
that a typical non-trivial periodic graph has intermedi- 
ate behavior with < p a < 1, as long as there is at least 
one gap and at least one band. One explanation of this 
phenomenon is that the graph of Fig. ITVa) (for example) 
can be viewed as a line with periodic (^-potential (the 
so-called Kronig-Penney model) whose strength is mo- 
mentum dependent. In such an analogue, the strength of 
the (^-potential oscillates between infinity and zero, which 
in effect alternates between disconnecting the graph and 
having a perfect transmission, resulting with the inter- 
mediate values < p a < 1. We refer the reader to [91125] 
for similar discussions. 



One can also consider dressing the network with a 
bounded periodic potential and/or changing the vertex 
conditions from the ones we considered. This should 
not affect our results qualitatively, as the influence of a 
potential or vertex conditions decreases in the k — > 00 
limit. However, this case is technically more difficult 
since the fc-dependence in equation (I5J) would become 
more involved. To overcome these difficulties, methods 
developed in |27l 128] might prove useful. 

Some further interesting spectral questions are now 
within reach. One may obtain bounds on possible sizes 
of bands (gaps) and deduce the specific edge lengths for 
which they are attained. Furthermore, the gap open- 
ing mechanism, a well studied subject on its own right 
[291 I3"D] , can be better understood by examining the sub- 
domains of the torus which do not intersect S. In ad- 
dition, the topological meaning of p a should be further 
investigated — does it relate to some other graph invari- 
ants or does it provide a brand new piece of information 
on the underlying graph? 

Finally, we make another step forward by extending 
the discussion to eigenfunction properties. The number 
of zeros of an eigenfunction was recently found to be con- 
nected with the stability of the corresponding eigenvalue 
with respect to magnetic perturbations [31 33 . The sta- 
bility is described by the Morse index of the eigenvalue 
and most strikingly, this Morse index can be shown to 
be a well defined function on the torus, not depending 
on the direction of the flow (i.e., on graph edge lengths) 
[34] , This leads to new and exciting findings on the dis- 
tribution of number of zeros of graph eigenfunctions |35| . 
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